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ABSTRACT 

O 

ON ■ We analyze the cut-off dependence of the fermion contribution to the finite tem- 

perature free energy density in 0(g 2 ) lattice perturbation theory for several improved 
staggered fermion actions. Cut-off effects are drastically reduced for the Naik action 
and an action with improved rotational symmetry of the quark propagator. We show 
qj | that improvement of rotational symmetry at 0(g 2 ) further reduces cut-off effects in 

thermodynamic observables. We also show that the introduction of fat-links does 
not have a significant influence on cut-off distortions at 0(g 2 ). 
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1 Introduction 



During recent years much progress has been made in the analysis of the thermo- 
dynamics of SU(N) gauge theories. Based on detailed studies of the cut-off de- 
pendence of thermodynamic observables in lattice regularized gauge theories an 
extrapolation of results for the equation of state to the continuum limit has been 
performed for the first time 0J. In addition, it could be shown that the influence 
of a finite lattice cut-off on thermodynamic observables is greatly reduced when im- 
proved actions such as 0(a 2 ) Symanzik improved or tree level perfect || gauge 
actions are used. By construction tree level improved actions do reduce the cut- 
off dependence of thermodynamic observables in the infinite temperature, ideal gas 
limit of gauge theories ||, They do, however, also seem to reduce cut-off effects 
at non-zero values of the gauge coupling g 2 , i.e. at finite temperature. Similar 
conclusions have been drawn from a first analysis of QCD thermodynamics with an 
improved staggered fermion action ||. 

Using tree level improved gauge and fermion actions it is now possible to perform 
lattice calculations at finite temperature, which at least in the infinite temperature 
limit lead to acceptably small systematic cut-off errors. Numerical simulations for 
SU (N) gauge theories as well as QCD with dynamical fermions || give some indi- 
cations that this improvement carries over also into the temperature regime close to 
the deconfinement phase transition. This suggests that improved actions also show 
a reduced cut-off dependence at non-zero gauge coupling g 2 . It would, of course, be 
nice to go on with the improvement program and construct systematically 0(g 2 a 2 ) 
improved or even 1-loop perfect actions which could be used for, e.g. thermody- 
namic, calculations. However, it is well known that this generates a large number 
of new terms, including 4-fermion operators, that would contribute to the 1-loop 
improved fermion action and would make such an action impractical for numerical 
calculations 0. We thus will aim at a more moderate goal in this paper. Within 
a class of tree level improved actions we will look for actions which lead to small 
cut-off errors in thermodynamic observables at 1-loop level. 

In the case of the standard staggered fermion action it is known that deviations 
from continuum perturbation theory are large for thermodynamic observables also 
at 0(g 2 ) ||. It is the purpose of this paper to analyze in 1-loop lattice perturbation 
theory in how far improved actions do reduce the cut-off dependence at this order. 
Besides analyzing tree level improved gauge and fermion actions we also discuss the 
influence of 1-loop improvement in the gauge sector on fermionic observables and 
construct a new fermion action with improved rotational symmetry at 0(g 2 ). Fur- 
thermore, we investigate in how far the use of fat links, which have been introduced 
to improve the flavour symmetry of the staggered fermion action ||, does influence 
the calculation of thermodynamic observables. Fat links modify cut-off dependent 
terms in thermodynamic observables at 0(g 2 a 2 ). 
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This paper is organized as follows. In the next section we introduce the various 
staggered fermion actions we are going to analyze in 1-loop perturbation theory and 
define the basic quantities needed in our perturbative calculations. In the following 
subsections we determine the coefficients of these actions to improve rotational sym- 
metry at tree-level and in 1-loop order. In section 3 we present results from a 1-loop 
calculation of the fermion contribution to the free energy density at finite tempera- 
ture. In section 4 we give our conclusions. Some further details of our calculations 
are summarized in the appendix. 



2 Improved fermion actions 

As starting point for our analysis we consider a generalized form of the naive fermion 
action consisting of terms which respect the hypercube structure of the staggered 
fermion formulation. In addition to the standard 1-link term this action also includes 
all possible 3-link terms resulting either from introduction of a higher difference 
scheme for the discretization of the fermion action or from the smearing of the 
1-link term (fat-links): 
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The fat-link improved one-link term, the linear 3-link term and the angular 3-link 
term appear with coefficients whose dependence on the gauge coupling we have 
parameterized as follows 
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Together with the fat-link-weight u there is thus a total number of 7 parameters in 
the action. In order to reproduce the correct naive continuum limit the coefficients 
have to satisfy two constraints 
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Our aim is to fix these coefficients such that the rotational symmetry of the fermion 
propagator is improved up to one loop order. 

In a first step we consider the free fermion propagator. Here we derive further con- 
straints for the tree level coefficients c\j. In particular we will consider the two 
cases where only one or the other of the three link terms contributes. Then the 
constraints fix the tree level coefficients to certain values. 

In the second step we calculate the self energy contributions to the fermion propaga- 

(2) 

tor. Here the one-loop coefficients c\j come into play and can be tuned to improve 
also the rotational symmetry in one-loop order. 

As gluon part of the action we choose the tree level improved lx2-action: 
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We note that also here the tree- level coefficients 5/3 and 1/6 could be further im- 
proved at 0(g 2 ). These corrections, however, are higher order corrections for the 
analysis of the fermion contributions to thermodynamic observables which we will 
present in the following. 

a Note that the N c dependence of the 1-loop coefficients has been factored out here explicitely. 
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2.1 Improvement of rotational symmetry at tree level 
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Figure 1: Dispersion relation E = E(p) with p = (p, 0, 0) for the standard staggered 
fermion action (blue), the Naik-action (green) and the p4-action (red) compared to 
the continuum dispersion relation E = p. Solid and dashed lines denote real and 
complex poles of the propagator respectively. 



The inverse free fermion propagator of the action |2.1| in momentum space takes 
the form 



A (or 1 / \ ~i E/x7/xfe/i(p) + m ~i Y.plv.Kip) + 
F E p h 2 p (p)+m 2 D(°)(p)+m 2 



m 



(2.6) 



where the functions h^(p) are given in appendix A.l| (Eq. |A.2|) . 
In order to obtain a rotational invariant free propagator we should achieve that 
D^\p) = J2u h fJi (p)h IJi (p) is a function of p 2 only. For the standard one-link prop- 
agator of the staggered fermion action this requirement for a rotational invariant 
propagator is violated at 0(p 4 ). Expanding D^(p), i.e. the trigonometric functions 
appearing in h^p) (see appendix |A.1| ), in orders of p a further constraint to the co- 
efficients Cj°i can be determined such that the free propagator is rotational invariant 



up to order p 4 , for details see appendix |A.1| . The constraint is 
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Two simple choices which eliminate the bended or straight three link terms in |2.1| , 
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Figure 2: Dispersion relation E = E(p) with p = (p/^/2,p/\^2, 0) for the standard 
staggered fermion action (green), the Naik-action (blue) and the p4-action (red) 
compared to the continuum dispersion relation E = p. Solid and dashed lines 
denote real and complex poles of the propagator respectively. 

respectively, are to set either cfl = which yields the familiar Naik-action 
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or to set C30 = which leads to 
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which we will call the p4-action. 

Since for the Naik-action the 0(p 4 ) terms are completely eliminated one gets even 
an 0(a 2 ) improvement, which seems quite accidental in this approach. 

A first impression of how the improvement works is obtained by inspecting the 
dispersion relation E = E(p x ,p y ,p z ) resulting from the poles of the free propagator 

for massless fermions , which are solutions of D(°\iE,p) = 0. Fi Kures ffl and 

show the dispersion relations for on-axis momenta p = (p, 0, 0) and momenta 
on the planar diagonal p = (p/y/2, p/y/2, 0) respectively for the standard fermion 
action, the Naik-action and the p4-action in comparison to the continuum relation 
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E(p) = p. In case of the Naik- action for on-axis momenta as in case of Naik- and p4- 
action for planar diagonal momenta there are also complex poles of the propagator. 
The real part of these poles is plotted as thin line for distinction to the real poles. 
In the continuum limit a — > only the branch in the lower left corner of the E-p 
plane survives as E = E cont ■ a and p = p cont ■ a. However, the finite cut-off a leads 
to deviations from the continuum dispersion relation on finite lattices. But the 
plots indicate that the dispersion relations for the improved actions are close to the 
continuum for a much larger range of momenta compared to the standard action. In 
particular the dispersion relation of the p4-action is close to the continuum relation 
in nearly half the Brillouin zone for on-axis momenta. As will be discussed in section 
|3| this leads to a strong reduction of the cut-off dependence in the ideal gas (g — > 0) 
limit. 



2.2 Improvement of rotational symmetry at 0(g 2 ) 

In order to reduce the violation of rotational symmetry also at 0(g 2 ) we again chose 
to look at the fermion propagator, as it is the most fundamental observable in the 
calculation of loop corrections to thermodynamic observables. 
Expanding the link variables U^(x) = exp(igA^(x)) one finds an expansion of the 
action in powers of g, 



S = + gS® + g 2 S^ + ^£(0) + Q ^ 



(2.10) 



Now expanding e s one gets the following contributions to the fermion propagator 
in momentum space up to 0(g 2 ) [] : 
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where the subscripts indicate that the expectation values are taken with respect to 
the free action and that only connected parts are taken into account. 



-5(0) 



(O} = Y o j [DU] [D^D^j] O e- s(0) , Z Q = J [DU] [D^Dtfj 

b Contributions arising from the expansion coefficients of the gluon action Sq" 1 and {Sq^) 
vanish because they lead to disconnected diagrams. 
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(4 1} ) 2 and are the self-energy terms arising from the parts independent of the 
1-loop coefficients c\j, whereas Sp^ denotes the part depending on these. We note 
that Sp^ and Sp^ are formally the same terms apart from an exchange of tree- 
level and 1-loop coefficients. The explicit form of these contributions are given in 
appendix |A.1| . 

An integration over the fermion and gluon fields yields, for massless fermions, 



^P'\p) 



E ^K 0) (p)+ff 2 (£ M (p)+D«(p))) 



A 



(o)- 



(p) 



with 
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+ ■ 



The explicit results are given in appendix |A.3| . The self-energy term contains, be- 



sides the usual continuum contribution (Sy) 2 , also a 1-loop term , Sp\ which is 
a pure lattice artifact. The term, which arises from the non-zero 1-loop coef- 

(2) 

ficients c\L acts like a counter-term in the sense that it ought to adjust the 1-loop 
contributions, although it does not remove any divergence. 

To achieve a rotational invariant fermion propagator to 0(g 2 ) we consider the de- 
nominator of A F (p), which can be written as 



V(p) = £ (D<M(p) + g 2 (E„(p) + D^{p)) 

= E^CP^CP) + 2 ^E4 0) (P) fa.(p) + Df\p)) + 0(g 4 



and demand that T> should be a function of p 2 only, analogous to the tree level 
case discussed in section EO. By this demand one obtains constraints for the 1-loop 

(2) 

coefficients c\ j analogous to the approach in the previous section. But in contrast 
to the tree level case a straight forward expansion in powers of p is not possible due 
to the logarithmic divergence in the contributions arising from (Sp^) 2 . 
Alternatively we consider on-axis and off-axis momenta having the same magnitude, 
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Figure 3: Solution of the equation |2.11| with complete contributions (blue) and 
without {Sp^) 2 contribution (red) for (a) the Naik-action and (b) the p4-action. 
Lines in case of the complete contributions are extrapolated by a polynomial fit. 



Pi = (ft, 0,0,0) and p 2 = (ft/y2,p/\/2, 0, 0), and then determine the coefficients by 
solving the equation | 



V(p x ) = V{p 2 



(2.11) 



for small momenta p. 

In general T>(p) depends also on the fat-link- weight u. But since the fat-link- weight 
has a negligible effect on thermodynamics, as it will turn out in the next section, 
where we look at the free energy density at 0(g 2 ), we fix the value of the fat-link- 
weight to u=0 for this calculation. The restriction to one or the other 3-link term 



in the action together with the constraint [2.4| leaves only one free 1-loop parameter 



which can be determined as solution of equation 2.11 



The term is a simple function of p. In the different parts of the Sp con- 
tribution the p dependence can be factored out of the gluonic integrals. Thus, they 



can be computed once for all values of p. In appendix A. 3 we show that only three 

(2) 

distinct gluonic integrals appear in the S F contributions. These gluon-integrals 
have been calculated to high accuracy using Gauss-Legendre integration. 
In principle one can expand these contributions in powers of p as in the previous 
section, but in contrast to that the (S^) 2 contribution can not be expanded due 
to a logarithmic divergence in the p — > limit. The gluon- loops have to be calcu- 
lated for each of the momenta p separately, since there is no factorization in gluonic 
and fermionic parts. The mentioned divergence causes also numerical difficulties 
for small fermion momenta p, where the pole of the inner gluon line, which is at 



The term ^u^^ip)^/! 1 (p) i s already rotational symmetric to order p A with the tree-level 



)(0) 



coefficients derived in section 2.1. Therefore only the term proportional to g 2 has to be taken into 
account. 
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p — k = 0, runs into the pole of the gluon propagator at k — 0. Therefore some care 
has to be taken to get reliable results. In practice this problem prohibits to perform 
our numerical integrations at arbitrary small p. 



Solutions of the equation 2.11 are plotted in Figure H for the Naik-action and 



(2) 

the p4-action as function of the momentum p. The squares are solutions for c\{, 
one obtains including only the term and the Sp^ part of S, i.e. setting (Sp^) 2 
to zero. Since the p dependence of these parts is known explicitely we get precise 
results down to very small momenta. The triangles represent the solutions including 
all contributions. The scattering of points at small p is due to the numerical difficul- 
ties in this regime for the (S { p ] ) 2 part. In the Naik case we obtain stable results for 
p > 0.1, whereas in the p4 case the results are trustable only for p > 0.3. In these 
regimes we have increased the number of Gaussian points until the values became 
stable, i.e. up to 320 Gaussian points in each direction for the small momenta. 
These plots indicate that the main contribution to the violation of rotational sym- 

(2) 

metry originates from the S F part of the self-energy term which is quite reasonable 
since this term is a pure lattice artifact. 

Using a polynomial fit we get the extrapolated values of the coefficients: 

Naik-action : c$ = -0.0085 , c$ = 0.00283 
p4-action : eg = 0.0165 , = -0.00275 

It is remarkable that these coefficients indicate an over- improvement using the tree- 
level improvement at finite g, as the tree- level coefficients are corrected towards the 
standard staggered action. 



3 Free energy in 1-loop perturbation theory 

Thermodynamic observables show large cut-off dependencies in finite temperature 
lattice simulations. A suitable quantity to study the influence of finite cut-off effects 
on thermodynamic observables is the perturbative high temperature limit which is 
known from continuum perturbation theory up to order g 2 for a long time ||10|| d . 



The free energy density for rif flavours of massless quarks and the SU(N C ) colour 
group is 



/cent = &, G + ftl, F + 9 2 (f ( 2t,G + fcl, F ) + 0(g 3 ) , where 
1 r 
45' 



d All perturbatively calculable terms (up to 0(g 5 )) have been calculated only recently for the 
pure gauge sector 
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/iol^/T 4 = n f ^-(Nl-l) . (3.1) 



Up to this order the thermodynamic relation e = 3p = —3/ holds, where e is the 

energy density, p is the pressure and / is the free energy density. 

On lattices with temporal extension N T large 0(1/N 2 ), i.e. 0((aT) 2 ), deviations 



from |3.1| are found in the gluonic contribution for the standard plaquette action as 
well as in the fermionic part for the standard Kogut-Susskind action. 
Tree level improved gauge actions reduce these cut-off effects in the leading order 
perturbative ideal gas limit. It has been shown that the improvement persists at 
finite temperatures well below those at which the ideal gas limit is approached 
0. In the fermionic sector thermodynamic observables show even larger deviations 
from the continuum Stefan-Boltzmann limit using the standard staggered fermion 
formulation ||. By construction the Naik-action also reduces these deviations in 
the fermionic contributions to the energy density e in leading order. 
At 1-loop order it is known that the standard staggered action also leads to large 
cut-off dependent deviations in the energy density e [|J. We will analyze in the 
following in how far these are reduced for the improved actions constructed in the 
previous section. 

We calculate the fermionic contribution to the free energy density fp in lattice 
perturbation theory up to order g 2 for various fermion actions introduced in the 
previous section including the 1-loop improved actions and fat-link improvement. 
Only at this order will the influence of fat-link improvement, 1-loop coefficients and 
also the chosen gauge-action show up. 

Since we are dealing with naive fermions the number of flavours is rif — 16. Up to 1- 
loop order lattice perturbation theory is equivalent for naive and staggered fermions 
and the correct nj dependence which only is a multiplicative factor at this order, 
can easily be introduced at the end [n|. We always chose the bare mass m = to 



get comparable results. 

The reason that we chose the free energy density / to look at is that it is defined 
by a very simple relation, just the logarithm of the partition function, instead of a 
derivative of the logarithm of the partition in case of the energy density e: 



v r 
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0(g 3 )\ .(3.2) 
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Here again we expanded the fermion action in powers of g. The subscript zero is 
defined as in the previous section and pure gluonic contributions are neglected also. 



Taking the logarithm of |3.2| one finds: 
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Explicit expressions for the expansion coefficients are given in appendix A.4 . The 
various contributions correspond to the following diagrams: 



f (o) 
If 




5(0) \ 

bp /o 




g i=> 




Obviously there is a correspondence between these contributions and those to the 
fermion propagator discussed in the previous section. One gets the graphs just by 
connecting the incoming and outgoing fermion lines of the diagrams of the fermion 
propagator. The Sp^ contribution includes the whole dependence on the 1-loop co- 
efficients and does not contribute if these are ignored. The Sp term again can be 
factorized in fermionic and gluonic loop integrals in contrast to the (Sp^) 2 term. 
In order to extract the finite temperature contribution we calculated the difference 
of the free energy density on N% x iV T and N* lattices. The spatial extension was 
set to infinity, N a = oo. The loop integrals over spatial momenta were done us- 
ing Gauss-Legendre integration. For the zero temperature contributions also the 
temporal momenta were integrated over with Gauss-Legendre integration, while for 
the finite temperature contributions the finite sums over temporal momenta were 
carried out explicitly. 

The numerical effort to calculate the 2-loop integrals in particular for the (Sp^) 2 
term is large and grows like (^Gauss-points) 8 , therefore the number of Gaussian 
points was at maximum 32 in each direction, which is a total number of about 10 12 
points. For small temporal extension N T < 8 one finds a good convergence of the 
2-loop contributions increasing the number of Gaussian points up to 32. For larger 
values of iV T , however, the limit is not reached with this maximum number of points. 
However, since the results show a clear (^Gauss-points) -4 behaviour they can be 
extrapolated to infinite number of Gaussian points also for iV T > 8. Of course the 
uncertainty of the extrapolation grows with increasing iV r , but since our main in- 
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Figure 4: Fermionic tree level contribution to the free energy density normalized to 
the fermionic continuum tree level contribution as function of temporal extension N T 
for standard Kogut-Susskind action (black), Naik-action (red) and p4-action (blue). 



terest lies in the deviations from continuum perturbation theory for small values of 
N T this does not affect our conclusions. 

Figure § shows the deviations of the fermionic tree-level contributions to the free 
energy density from the Stefan-Boltzmann limit for the standard staggered fermion 
action, the Naik-action and the p4-action. The improved actions reduce the large 
cut-off dependent deviation of the standard action drastically from more than 20 
at N T =A to about 20% in case of the Naik-action and about 40% in case of the 
p4- action. At N T =6 and greater the improvement works even better. Here the de- 
viations for the p4-action are below a few percent compared to nearly 80% for the 
standard action at iV T =6. 

The results for the fermionic part of the free energy density at order g 2 , f^ 2 \ nor- 
malized to the continuum order g 2 contribution, f^l^, are plotted in Figures [5] and 
|6|. Figure [| shows the influence of fat-link improvement and the chosen gauge action 
for the considered tree-level improved and the standard staggered fermion actions, 
i.e. 1-loop coefficients are set to zero. For all actions large deviations from the con- 
tinuum value are found at small values of N T , although the p4-actions approach the 
continuum value most rapidly. The main differences are caused by the choice of tree- 
level coefficients of the fermion action rather than the gauge action or the fat-link 
parameter This is not unexpected, because fat-link improvement affects the in- 

°The values of u> are chosen in the same range as used by the MILC collaboration for investi- 



13 




Figure 5: Fermionic 1-loop contribution to the free energy density normalized to 
the fermionic continuum 1-loop contribution as function of temporal extension N T 
for various fermion actions. Different colours represent the value of the fat-link 
parameter uj = 0.0, 0.2, 0.4. Circles denote the use of the standard gluon action, 
open squares the use of the 1x2 gluon action. In the Naik- and p4-action only tree- 
level coefficients are included, i.e. cfj = 0. The values for N T > 8 are extrapolated. 

frared behaviour and reduces flavour symmetry breaking rather than the ultraviolet 
regime which is responsible for the high temperature behaviour. Further, improved 
gauge actions are designed to reduce cut-off distortions in the gluonic contributions, 
but they have shown to reduce flavour symmetry breaking also , which indicates 
rather an influence on infrared sensitive fermionic observables. 

The influence of 1-loop improvement of Naik- and p4-action provided in the 
previous section on the high temperature behaviour at order g 2 is shown in Figure 
Here in all cases the 1x2 gauge action is used and the fat link parameter equals 
zero. The plot indicates a large reduction of cut-off effects at order g 2 for both types 
of 1-loop improved fermion actions compared to the corresponding tree-level action. 
Already at iV T =6 the deviations are reduced to less than 50% for the 1-loop p4- 
action and to about 140% for the 1-loop Naik-action, where the 1-loop contribution 
for the standard staggered action is about 4 times larger than the continuum 1-loop 
value. 

As stated in the previous section, the 1-loop coefficients correct the tree-level 
coefficients towards the standard action in both cases, Naik- and p4-action. In 

gation of flavour symmetry breaking Q . 
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Figure 6: Fermionic 1-loop contribution to the free energy density normalized to the 
fermionic continuum tree level contribution as function of temporal extension N T for 
the Naik-action (red) and the p4-action (blue) including only tree-level improvement, 
i.e. c[ 2 j = 0, (open circles) and including also 1-loop improvement (squares). Always 
the 1x2 gluon action is used. The values for N T > 8 are extrapolated. 

(2) 

correspondence to that we observe a change of the sign of f F comparing standard 
and tree-level improved actions, i.e. an "over- improvement" concerning the order g 2 , 
which is reduced when including the 1-loop coefficients. 



4 Conclusion 



In this work we have presented a perturbative construction of staggered fermion 
actions with improved rotational symmetry of the fermion propagator at tree-level 
as well as in 1-loop order. In particular we considered tree-level and 1-loop improved 
versions of the Naik- and the p4-action both with and without additional fat-link 
improvement. 

The tree-level improvement shows up in the dispersion relation which is close to the 
continuum relation in a much larger momentum range for the improved actions than 
for the standard action. 

A perturbative calculation of the fermionic tree-level and 1-loop contributions to the 
free energy density / has been performed to investigate the influence of improved 
rotational symmetry and fat-link improvement on the high temperature behaviour. 
At tree level as well as in 1-loop order the standard staggered action leads to large 
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cut-off dependent deviations from the perturbative continuum value. We find a large 
reduction of these ultraviolet cut-off effects at tree-level for the tree-level improved 
actions. For example at N T = 4 the remaining cut-off effects are about a factor 10 
smaller than the cut-off effects of the standard action. 

Of course the coefficients of the 0(a 2 ) deviations depend on the observable. An 
earlier calculation of fermionic contributions to the energy density e at tree-level 
has shown a reduction of cut-off distortions for the p4- and the Naik-action which is 
comparable to the results for the free energy density we present in this work although 
in case of the energy density one finds even smaller deviations for the p4-action of 
about 1% at N T =4 §. 

Tree- level improvement, however, does not help at 1-loop order as expected. 
Our analysis showed that cut-off distortions are of comparable magnitude for the 
standard action as for the tree-level improved actions in the 1-loop contributions to 
the free energy density. It also turned out that fat-link improvement has a rather 
negligible influence on the high temperature behaviour, i.e. it neither reduces nor 
increases the cut-off effects significantly. This is not surprising since fat-link im- 
provement was designed to affect the infrared rather than the ultraviolet behaviour. 
Finally we find a large reduction of these ultraviolet cut-off effects at 1-loop order for 
the actions with improved rotational symmetry up to 0(g 2 ). Already at iV T =6 the 
deviations are reduced to less than half of the continuum 1-loop value for the 1-loop 
p4- action and to about 1.4-times the continuum value for the 1-loop Naik-action, 
where the 1-loop contribution for the standard staggered action is about 4 times 
larger than the continuum 1-loop value. 

Therefore the 1-loop improved p4-action in combination with fat-links is a particu- 
larly suitable candidate for finite temperature simulations in full QCD, combining 
the advantages of an improved high temperature behaviour up to 1-loop order with 
an improved flavour symmetry. And last but not least it retains the simplicity of 
the contributing local operators. 

Acknowledgments This work was partly supported by the TMR net- 
work Finite Temperature Phase Transitions in Particle Physics, EU contract no. 
ERBFMRX-CT97-0122. The work of U.M.H. was partly supported by DOE con- 
tracts DE-FG05-85ER250000 and DE-FG05-96ER40979, and he gratefully acknowl- 
edges support from the Zentrum fur Interdisziplinare Forschung, Universitat Biele- 
feld. 



16 



A Appendix 



We give here explicit expressions as well as some technical details of the perturbative 
calculations. It is organized as follows: In part [A.l| we define the fundamental terms 
and functions we refer to in the following parts. Part A.2j and part A. 3 contain some 
details of the determination of the tree-level and 1-loop coefficients respectively. 
In part A. 4 we give the explicit results of the 1-loop calculation of the fermionic 
contribution to the free energy density. 



A.l General 1-loop results and definitions 

For the trigonometric functions we use the short hand notation: 



s M (p) = sra(p„) 
ty(p) = cos(p M ) 



The inverse gluon propagator of the 1x2 action was calculated in ||. For the 
standard plaquette action as for the 1x2 action the propagator can be written as 



Ag^O) = D Gli {k)8^ v - E GlhV {k) + f g^k)g u (k) , 
D Gfl (k) = 4 01,1 £ 4(k/2) - 16 oi, 2 £ 4(k/2) (2 - £(k/2) - s%k/2)) , 

V V 

E G „A k ) = 4 s»{k/2)s v {k/2) 

-16 a 1>2 s^k/2)s v {kl2) (2 - s 2 v {k/2) - s%k/2)) , 

with ai i = 1) a i.2 = for the standard Wilson one-plaquette action and a^i = 5/3, 
ai ; 2 = —1/12 for the 1x2 action. In both cases we choose the gauge fixing term 
9fj,(k) = 2s At (/c/2) and Feynman gauge ^ = 1. 

The basis of our 1-loop calculations is an expansion of the fermion action in 
powers of the bare coupling g: 

S F = SP + gS§> + g 2 sf + g 2 §P + 0(g 3 ) . (A.l) 

To achieve that, we expand the exponential representation of the link variables 
U^(x) = exp(igaA^(x)), with the gauge fields = Yl^=i 1 A b ^ b and normalization 
2Tr A a A 6 = 5 a b for the group generators A a . The lattice spacing is set to a = 1. In 
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momentum space we find then 



q(0) 



s 



(0) 



7(1) 



s 



(2) 



'^(p) \J2^K(p) + m \ ^(p) 



/V 2 — 1 /• 



4/7 / $(P) E7^f ) (P,*i ) M#-fei-fc2 



where 



J2 K KP (p, k) A p {k) , 
p 

^ L^ p , a (p,ki, k 2 ) A p (k l )A a (k 2 ) 



p,cr 



Here we use the following definitions also referred to in appendices |A.2| - |A.4j: 



Mp) 
^ 2) (p) 

L p ;p,a(P' k l,k 2 ) 



2 s» [eg + 2 eg C -( 2 P)J + 2 4>m(3p) , 
2 s» [eg + 2 eg J] C -( 2 P)] + 2 4 2 k(3p) , 



eg -4J>;P, fc) + 45 A^\p, k) + eg k) , 



(0) 



p;p 



(A.2) 



eg P, fci, h) + eg Bg£ (p, k h k 2 ) + eg Bg# (p, fei, 



2J , 



where 



= 2c>- k/2) 

<W - (<W £ s 2 u (k/2) - (1 - s „(*/2) a ,(*/2) 

= 2^-^/2-^/2) 

<W<W - ; ( <W<W s U k i/2 + k 2 /2) 

1 + bu; V ^ 

+2 (1 - <^ p ) <f„ i(y s M (A; 1 /2)^(A; 2 /2)c p (A; 1 /2 + fc 2 /2) 
-<W (1 - <W) K(V 2 ) - ic p {k 2 /2)} s a (ki + k 2 /2) 

-<W (1 - 6 P}P ) {s p (h/2) + ic p (h/2)} s a (k 2 + h/2) 
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-- 2^, pC/i (3p-3fc/2)(3-4sJ(fc/2) 
= 2 5^ 6^ s p (3p - 3h/2 - 3k 2 /2) [9 - 4sJ(fci/2 + k 2 /2) 

-4^(^/2) - 4 S J(A; 2 /2) - 4^(^/2 - fc 2 /2) 
= 4 C/i (p - k/2) <5 M , P X! - fc)c^(A;) 

-8 (1 - <^ p ) s> - A;/2)c M (A;/2) Sp (2p - k)c p (k/2) , 
-- 4 — ki/2 — k 2 /2) 

<W c v(k\ + k 2 )c u (2p -h- k 2 ) 

+2 (1-^)^^/2 + ^/2) 

■{c p (2p - 3h/2 - k 2 /2) + c p (2p - fci - fc 2 )c p (A;i/2 + fc 2 /2)) 

+8 (1 - <W) c pG° - fci/2)c T (A;2/2)s ff (2p - 2 A* - fc 2 ) 
+8 V (1 - <WKGo - fci - fc 2 /2)c p (A; 1 /2)s p (2p - £4) . 



The integral symbols with subscripts for fermion and gluon momenta, denoted 
by the letters p and k respectively, actually denote a sum over the momentum modes 
on finite iVj x N T . They are defined as: 



f wit J k, = 2 -^ , -\N a <n p <\N a -l, 

J* N%N T ^ \ k A = ^ , -\N T <n 4 <\N T -l 

f 1 ^ ..,/ = ^ , 0<n M <7V CT -l, ^^4 

A " 1 P4 = ^, 0<n 4 <iV T -l 



The shift of temporal fermion momentum modes, p 4 , represents antiperiodic bound- 
ary conditions in the time direction. As we consider the infinite volume limit the sum 
becomes an integral over the Brillouin zone, l/N a J2ii ~^ l/(2vr) f^d/j, % = 1,2,3, 
for both fermionic and gluonic momenta, / = p, k. The same holds for the time 
direction in the zero temperature contributions. 

The path integrals over the fermionic and gluonic fields are evaluated using the 
identities 

Al(k)Al(-l)) o = ^1 A ^„(fc) S k , 5«' b , 

r*tp)m)) Q = - — ^ — — ^ ^ > 
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with 



m 2 



A. 2 Tree level coefficients for improved rotational symme- 
try 

An expansion of the free inverse fermion propagator up to order p yields 

d (0 \p) = E = E H A [ A + 2B ^l + 2 ^ E Pt ) + °(p 6 ) . 



with coefficients 



A = 2cg> + 12cg + , 
St - -ic (0) - 2c (0) - 9c (0) 

5 2 = -*n . 



Obviously the condition B\ = B 2 has to be satisfied to achieve rotational symmetry 
up to order p 4 , which leads to the constraint 

JO) , 2? (0) fi (0) _ _ (0) 
c 10 -\-zi c 30 -f- o c l2 — &i c 1)2 • 

The expansion shows that setting c[°2 = 0, i.e. the Naik- action, leads to vanishing 
0{p A ) coefficients B± — B 2 = 0, which corresponds even to an 0(a 2 ) improvement. 



A. 3 Improved rotational symmetry at 0(g 2 ) 

The explicite expressions for the 0(g 2 ) contributions to the fermion propagator are: 

D^(P) = -Kip) , 
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N 2 — 1 r 1 r 

Mp) = C 2N J k - ( k) E E [Mp - k)K v . p (p, k)K v .- p {p - k, -k) 



2N C Jks F (p-k)^r Ptp 

-2h u (p - k)K KP {p, k)K v , p {p - k, -k^A G \ rp (k) , 

A'- - I 



N 2 — 1 r 

= —fif- J k E L ^-M K -k)A G \^(k) . 



(2) 

For Cfj,(p) which originates from the S F term the p dependence can be factored out. 
To achieve this one has to express the trigonometrical functions of sums of p and 
k as sums of factorized ones. Then taking into account that A^ 1 (k) is an even 
function of all /c-directions and that A^ 1 (k) is an odd function of k^ and k u and 
an even function of the remaining /c-directions one finds that a lot of terms are odd 
functions in some /c-direction and therefore vanish upon integration. The remaining 
terms, without using fat links, i.e. taking oj = 0, include only three gluonic integrals: 



= -^^{s,{p)cfl2X l + s^p)c'il2X 2 



+ s,(p)J2cu(2p) cS J(4 J 1 + 16 J 3 )| , 



where the gluonic integrals are defined as 

J 2 = j k ^-%s\{k/2)-As\{k))A- G \^k) , 

J 3 = /(^(^A^^-^^c^^^^A^ 1 ^^)) . 

A calculation of these integrals to high accuracy using Gauss-Legendre integration 
gives 

X l = 0.1282908 , 
J 2 = 0.0752730 , 
J 3 = 0.5074740 . 
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A. 4 Contributions to free energy density 



The explicite expressions for the fermionic contributions to the free energy density 
up to 0(g 2 ) are: 



/» 



s 



(o)\ 
F /o 



-2N C I In [s F (p)] , 
Jp 

2{N 2 C - 1) 



E 



P Jk s F (p)s F (p-k) ^ K 
<W E h p(p) h p(P -k)+m 2 ) - 2h tl {p)h p {p - k) 



1 

P JkS F (p)^ jKj , 

(iv c 2 -i) [j2 h ^ p)h( " ){p) 



sP) = {Ni -i) J J J—y, Kip) E Wp> *, 



sj?(p) 



In (S^y the integration/sum can again be factored, similarly as for C^ip) in 



o. 
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